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1 Introduction 



The Gaudin model is an example of an integrable system in statistical quantum mechanics, as- 
sociated to N spin particles with magnetic interaction, based on N independent and commuting 
Hamiltonians. As originally formulated, this model is related to the Lie algebra SL{2). In this 
simplest case, the problem of the diagonalization of the Hamiltonians was partially solved. Us- 
ing a constructive method, known as Bethc Ansatz, Gaudin |^ constructed a family of common 
eigenvectors for all Hamiltonians. 

More generally, integrable systems can be associated to any semi-simple complex Lie algebra. 
For such systems, generalizations of the set of Gaudin Hamiltonians have been constructed.]^ The 
methods proposed for diagonalization of these Hamiltonians revealed remarkable connections be- 
tween integrable models and two dimensional conformal field theories. In this approach the 
structure of the Bethe vectors for the Gaudin model is related to the representation theory of affine 
Lie algebras. The diagonalization of Gaudin Hamiltonians is based on a concept of invariant func- 
tionals (correlation functions) on tensor products of representations of an affine Lie algebra at the 
critical level. 

Another approach for the problem of diagonalization, || is related to the connection be- 
tween eigenvectors of the Gaudin Hamiltonians and the solutions |^ of the Knizhnik-Zamolodchikov 
equations. In this approach, the common eigenvectors are constructed inductively and each eigen- 
vector leads to an integral solution of the KZ equations. Conversely, it was proved § that in the 
quasi-classical limit, the first term of the asymptotic solutions of KZ equation leads to a common 
eigenvector of the Hamiltonians. In this case the Bethe equations for the Gaudin model appear as 
conditions of critical points. 

In section 2 we review some aspects concerning the Gaudin model, originally formulated for the 
Lie algebra SL{2). We recall the expression of the N Hamiltonians, introduced by Gaudin. The 
space rj of physical states is defined as a tensor product of N finite-dimensional highest weight 
representations of the Lie algebra SL{2). It is decomposed as a direct sum of weight subspaces, i.e., 
O = ©Vm- We explain how the Bethe Ansatz works to construct common eigenvectors for Gaudin 
Hamiltonians in each eigenspace Vm- However, this method cannot give all the common eigenvectors 
of Gaudin Hamiltonians. Therefore, supplementary common eigenvectors have to be determined. 

In section 3 we give a general, recursive method to construct a basis of common eigenvectors in 
each invariant subspace Vm- Knowing a basis of Vm-i, we construct a family of linear independent 
common eigenvectors in Vm, which are nonsingular. This family is completed to a basis of Vm, 
by a basis of the subspace of singular vectors of Vm- Using an analog of the Gordan operator, 
we give a basis of Sing Vm- Finally, we show that the Bethe Ansatz gives a family of singular 
common eigenvectors for Gaudin Hamiltonians, but in some particular cases, this family could not 
be complete. 

Section 4 is devoted to a discussion of the generalized model associated to an arbitrary simple 
Lie algebra. We recall the generalization of the Gaudin Hamiltonians and explain how the Bethe 
Ansatz was generalized to construct common eigenvectors for these Hamiltonians. We show that 
Bethe equations, which appear as conditions that Bethe type vectors be common eigenvectors, are 
also conditions of singularity for them. This constructive method does not ensure the completeness 
of the system of eigenvectors. 

2 Gaudin model and the SL{2) Lie algebra 

The Gaudin spin model is related to the Lie algebra SL{2), with generators E, F and H satisfying 
the commutation relations: 

[E, F] = H; [H, E] = 2E; [H, F] = -2F 

For this Lie algebra we consider N finite-dimensional highest weight modules: Vx^^,..., Vx^ with 
highest weights Ai,..., A^v and highest weight vectors v\-^,..., vx,^. The tensor product of these N 
modules constitutes the space of physical states for a system of N spin particles: 

n^Vx,(E) .■.(E>Vx„. 



1 



The Lie algebra action on each vector v = vi ® ... 'E) vn of this tensor module is defined as 



N 



(1) 



where X^^^ denotes the operator on O which acts as X on the module and as the identity operator 
on all other factors: 



(g) ... Wat = wi 



)Xvi 



'Vn, 



(2) 



For such a system of N spin particles, Gaudin proposed a set of A'' Hamiltonians defined on fl, 
depending on A'^ distinct, complex parameters zi, ...,zn ■ 



^ 1 

ni{zi,...,ZN) = Yl Zi-Z- 
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V i = 1 N. 



(3) 



All these operators commute: 



[n,,Hj]=0, V^,J = l,...,iV, 



but they are not independent, because J2f=i = 0- can easily be verified that among the N 
Gaudin Hamiltonians there are exactly N — 1 which are independent. To integrate a system of N 

spin particles, with N degrees of freedom, the family of commuting operators is completed by the 
Cartan generator (also called the total spin operator), of which action on the tensor module is 



JV 



and which commutes with all Gaudin Hamiltonians. 

For this family of N independent Hamiltonians which commute there is a complete system of 
common eigenvectors in fl. Our goal is to construct a basis in Cl, of common eigenvectors of Gaudin 
Hamiltonians. 



2.1 The structure of the space 

In order to analyze the structure of f2 it is useful to recall some elements of the theory of highest 
weight representations of the Lie algebra SL{2). Such a representation is completely determined by 
a highest weight vector vx, on which the action of the algebra is given by 

Hvx = Xvx, Evx = 0, 
and the representation space is generated by vectors: 

{Vn = WAlneN • 

The action of the Lie algebra on these vectors is 

Hvn = {X-2n)vn, (4) 
Evn = n{\-n+l)vn-i, (5) 

FVn = Vn+l- (6) 

If the highest weight A is not a positive integer, the representation is infinite dimensional and 
irreducible. If A G N, then in {w„ = F" WAjnGN there is an invariant subspacc, generated by 
{ly'A+i, v\+2, ■■■} and the quotient representation is irreducible and finite dimensional, of dimension 
A + 1, generated by vectors {vn = -F" vx}n=o,...,\- We denote Vx this quotient representation, for 
which we have F" = for all n > A + 1. 

We will consider the space n as a tensor product of finite-dimensional representations Vx, which 
is completely determined by the vector vq = vxi (g) ... (8 vxf^, called vacuum vector and is generated 
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by vectors {F"^uai F^'^v\^... (g) F"^" v\^}ni=o.,...,\i- i=i....,N, with rii operators F apphed on the 
i*^ component. Such a vector can be written as a product of m — ni + n2 + ■■■ + un operators of 
type (||), denoted 

with 1 < fci < ... < kjn < N. Note that for finite-dimensional representations, m can vary between 
and a maximal value m^^^ = Si!=i ^i- 

Note also the particular action of the Hamiltonians on the vacuum vector: 

XiXj 



/ AT N / N 

^^0=i^Xk\vo, 'H^Vo=\- ^ 



(7) 



The Cartan operator H has a privileged place in the family of commuting Hamiltonians. Ac- 
cording to it gives a gradation of the representation spaces V\^, which induces a gradation of 
the tensor product module, on H-invariant subspaces: 



m—0 



where Vm is a weight subspace, of weight ^ 2m, also called space of spin deviation m. It is 

generated by H-eigenvectors: 

{vi^;^'-''^^) = F^'''l..F^''^-^vo}i<k,<...<k^<N, (8) 
which we call states of spin deviation m: 



'Wo- 



Note that for spin deviations m < min{Ai, Ajv}, a subspace Vm has the dimension 

dimV,„= E 1= E l = ^]v"+™-i- (9) 

l<fei <...</£„, <Af l<ji<...<j,„<Ar+m-l 

From the explicit form of the Gaudin Hamiltonians (^) and from the action (^)-@ of the SL(2) 
Lie algebra on the spaces , it follows that the weight subspaces Vm are invariant under the action 
of any Jii. Therefore, in each subspace Vm, we can construct a basis of common eigenvectors of Tii. 

2.2 The construction of common eigenvectors by Bethe Ansatz 

Bethe Ansatz is a method to construct a family of common eigenvectors for Gaudin Hamiltonians 
in each invariant subspace Vm , but this family is not a basis of Vm . 

Since Vm is generated by (0), any common eigenvector in Vm has the form 



N N 



kl=l fcm = l 

The central idea of the Bethe method is to consider the coefficients Ck-^,,,k^ as rational complex 
functions: 

1 1 
Cki...k„^ = ■•■ , 

Wl - Zfej Wm - Zk^ 

depending on some unknown, distinct, complex parameters: wi, ...,Wm- We call such a vector a 
Bethe vector. Hereafter we shortly present this method. 



3 



2.2.1 Bethe vectors of spin deviation m — 1 

The eigenspace Vi is generated by vectors {-^'^''^^0}^.^^ jy. A Bethe vector of spin deviation 
m = 1 is defined as an expansion on this basis, with rational coefficients depending on one complex 
parameter w: 

where we denoted by J-^{w) the operator on fl: 

N 

' W - Zk 



^ 1 

= Er73^^^'^- (10) 



k=l 

Straightforward calculations give the commutator 

[w-z^) {w-Zi)^^W-Zk 

Applying this operator on vq and using we obtain the action of a Gaudin Hamiltonian Hi on 
i>i{w): 



(l ^ 



^iM" E 



-Wo, 



Zi- Zj W - Zi j \^ W - Zk W - Zi 

j \k—l ' 

and we have the following lemma. 

Lemma 1 Given N distinct complex numbers {zi}^^^ jy and fixed, positive, integer highest weights 
N' Bethe vector ipiiw) of spin deviation m = 1 is an eigenvector for all Gaudin Hamil- 

tonians: 

Hiipiiw) ^ si -ipiiw), Vi = l,...,7V, 
if the complex parameter w satisfies the condition 

E - 0' (12) 

- Zk 

called Bethe equation associated to Vi. The eigenvalue sj of Hi, depends on the solution w of this 
equation: 



E 



\i \j Xi 



2 ^ , Zi ~ Zj w — Zi 
2.2.2 Bethe vectors of spin deviation m — 2 

The eigenspace V2 is generated by vectors {F'"''^^ F'^'^^^vq} ^^^^^j^^^j^- A Bethe vector of spin devi- 
ation m = 2 is defined as an expansion on this basis, with rational coefficients depending on two 
complex parameters wi, W2 ■ 

i>2iwi,W2) = T{wi)T{w2)vo. 

Remark that the order of the two operators J- is not significant because they commute. 

In order to compute the action of a Hamiltonian on this state, we have to commute 7i,; successively 
with the two operators J-. Using (|ll|) we obtain the commutator formula: 



2 / p(0 

VUi - W2 \Wi - Zi W2- Zi ' 
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From (|11|) and (y_3|) the action of a Gaudin Hamiltonian Tii on the state ■02 is 

'Hitp2{wi,W2) = I ^ ^ ^! ^ ^ \ 1p2iwi,W2) 

\ 2 . 'T-' . - Zj wi - Zi W2- z, 

- V + ^I^1{W2) 

Vfell '"'2 ~ •^'^ Wl - W2 J W2 - Zi 

and we have the following lemma. 

Lemma 2 T/ie Bethe vector i'2iwi,W2) of spin deviation m = 2 is an eigenvector for all Gaudin 
Hamiltonians: 

Hi'lp2iwi,W2) = sf V'2(wi, W2), V i = 1, ...,N, 
if the complex parameters wi , W2 satisfy the Bethe equations associated to V2 • 

y^^ + ^_^0, V^^ + ^^=0. (14) 

j-^ Wi~ Zk W2 - Wl i-^W2~ Zk Wi - W2 

The eigenvalues sf depend on the solutions wi,W2 of these equations: 

N 



XiXj A j 



2 . . Zi - Zj Wi - Zi W2- Zi 

2.2.3 Bethe vectors of spin deviation m 

In the subspace Vm, generated by vectors {F^''^\..F^'''^^vo} ^^1^ ^ a Bethe vector of spin de- 

viation m is defined as an expansion with coefficients depending on m complex parameters wi, Wm '■ 

1pm{wi, ...,W.m)vo = T{wi)...T{w„i)vQ. (15) 

The action of a Gaudin Hamiltonian Tii on this state is calculated by induction on m: 

Hi'lpm{wi,...,Wm)vQ = HtT{Wm)lJjm-liwi,...,Wm~l)vo 

= T{Wm)'Hi1p,n-l{wi, m„i_i)wo + i^m-i{wi, m„i_i) [Hi , T{w 

m )\ ^0 

+ [[Hi , T{Wm)] ,1pm-liwi, W„i_i)] Vo. 



(16) 



The first term is computed from the induction hypothesis, the second from (|11|) and the last term 
from relation ([l3|). Putting all together it follows that 

/ 1 ^ AiA. ™ A, \ 
J^/^ A, ™ 2 \ , , „ , 

-V V + V Vm-l(--,Wfc,--)«0, 

,tt - iJ^^k ^i-^kj Wk- Z, 

where ipm-ii-'-^Wk, ■■■) denotes (wi, Wfc+i, Wm). Hence, we obtain the following 

theorem. 

Theorem 3 The Bethe vector ^j^niwi, ...jWm^VQ of spin deviation m is a common eigenvector for 
all Gaudin Hamiltonians: 

Hi 1pm{wi, ...,Wm)vo = S™ V'mCwi, ...,Wm)vo, V i = 1, 
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if the complex parameters wi, Wm satisfy the Bethe equations associated to Vm-' 

AT m 

E^^+ E ^^=0, Vfc = l,...,m. (17) 

The eigenvalues s™ depend on the solution of these equations: 

N 



As observed in Ref . , the Bethe Ansatz does not give all the common eigenvectors of Gaudin 
Hamiltonians. Therefore, supplementary common eigenvectors have to be determined, which will 
be done in the following section. 

3 A construction of a basis of common eigenvectors for Gaudin 
Hamiltonians 

In this section we give a general method to construct a basis of common eigenvectors in each 
invariant subspace Vm- This is a recursive method, which will be applied for subspaces Vm, with 
m < min{Ai, Aat}. Knowing a basis of Vm-i, we construct a family of linear independent common 
eigenvectors in Vm, which are nonsingular. This family is completed to a basis of Vm by a basis of 
the subspace of singular vectors of Vm • 

There are two important properties of Gaudin Hamiltonians which are useful in this section, 
namely they commute with the operators E and F, on the tensor product module: 

[H„E] = 0, (18) 
[7^„F] = 0, (19) 



where E and F are defined by ([|) as 



N N 

E = ^£;«, F = ^F«. 



Definition 4 We call v" E fl a singular vector of if the generator E acts trivially on w^, 

Ev" = 0. 

We denote Sing the subspace of singular vectors in V. 

It was shown in Ref. that the dimension of the subspace Sing Vm is 

dim (Sing Vm)=C;;^+^_2. (20) 

Since according to (^), for m < min{Ai, Xn} the dimension of the space Vm is C™^jy-i, it follows 
that a basis of Sing Vm can be completed to a basis of Vm, by a family of C™^^^_2 nonsingular linear 
independent vectors. The space spanned by this family of nonsingular vectors is denoted NonSing 
Vm- Hence, 

Vm = Sing Vm © NonSing Vm- 
From the property ( |l8| ) we obtain the following. 
Lemma 5 Sing Vm is a C^^_^j^ _2- dimensional vector subspace ofVm o,nd it is Tii -invariant: 

(SingVm) C SingVm, = 1, iV. 
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Then Sing Vm admits a basis formed by singular common eigenvectors of Gaudin Hamiltonians TLi. 
Denote Bi^ this basis. 

From the property ( [l9| ) of the Gaudin Hamiltonians we can construct recursively a basis of 
common eigenvectors in NonSing Vm'- 

Proposition 6 The basis of common eigenvectors in NonSing Vm is obtained by the application of 
the operator F on all common eigenvectors (singular and nonsingular) which form the basis of the 
invariant subspace Vm-i- Denote B^^^ this basis: 

B'^^^YiBm-i). 
where Bm~i is a basis of common eigenvectors ofVm-i- 

Proof: (1) The vectors of B^ are linear independent. Indeed, for m < min{Ai, Ajv} the 
operator F : Vm-i — > Vm is injective, since KerF — {0}. Moreover, the injectivity of F implies 
that the number of elements of the family B^ is equal to the dimension of Vm-i, namely C™^^_2- 

(2) The vectors of B^ are common eigenvectors of Hi. If Vm-i G Sm-i is a common eigenvector: 

HiVm—l = O^m— 1^™ — li 



then, from the property (19) it follows that the vector Vm — Fwm-i G B^ is also an eigenvector, 
with the same eigenvalue am_i. 

(3) The vectors of B^ are nonsingular. 
Note first that in Vo there is only one vector, the vacuum vector vq. Thus, there is one common 
eigenvector vi =Fvq in the family S"'*, which is indeed nonsingular: 



Evi = EFwo = Hvo 



/ N \ 
\i=l ) 



Now suppose by induction that for all fc < m — 1, with m < min{Ai, A^v}, we have constructed 
a basis of V^ of the form U B"^^ , where S^* is obtained by application of the operator F on vectors 
of the basis Bk-\ of Vfe_i. Then for a vector of S"* we have 



N 



E = E F Vra-i = H Vm-i + F E Wm-1 = (^51 ~ " 1) j + ^ ^ '^'"-1- (21) 

If Wm-1 is singular, then the last term of ( pl| ) vanishes and v'^-^ is nonsingular, since 
^A, -2(to-1) >iVTO-2(m-l) >0, VA^ > 2. 

If Wm-i is nonsingular, then by the induction hypothesis, there is a fc G {1, ...,m — 1} such that 
Vm-\ = F''^m-i-fc '^ith ^m-i--fc ^ commou singular eigenvector of i3m-i-fc- Then relation ( |2l] ) 
becomes 



= (fc + 1) A, - 2(m - 1) + fc j Wm-i ^ 0, ViV>2. (22) 

4) The vectors of are linear independent of the singular vectors of the basis B^- 
Consider a null linear combination of the vectors v^ii) £ B^ and vf^{j) G B^'- 

^a(z)<f(z)+^6(j>^(j-) = 0. 

* 3 

Applying the operator E, the second sum (of singular vectors) vanishes. Hence, 

^a(^)E<^(^)=0. 
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It follows from ( p2|) that each vector Ei;"*(i) is colinear with the vector which form the 

basis of Vm-i- Then, the coefficients a{i) are all zero. Since v'^{j) form also a basis, in Sing Vm, 
the coefficients b{j) vanish too. 

Properties (l)-(4) from above show that B""^ — F {Bm^i) is a basis of NonSing Vm and then 
B"^ U B!^ is a basis of common eigenvectors in Vm- In this basis, the sub- family is completely 
determined by Proposition ^ whereas for the subset of singular vectors Bf^ we have only the existence 
Lemma |^ and not the structure. In the next section we will characterize the singular vectors of a 
tensor product module of SL{2). 

3.1 Singular vectors of the SL(2) tensor modules 

For a finite-dimensional highest weight module V\ of the Lie algebra SL{2) there is only one singular 

N 

vector: the highest weight vector v\. Nevertheless, for tensor product modules = ®V\- the 

1=1 

N _ 

vacuum vector vq = ® vx is not the only singular vector. There exist singular vectors in every 
subspace Vm- 

3.1.1 ThecaseA^ = 2 

Consider first N = 2 and Q = Vx^ <SiV\2 with the decomposition on invariant subspaces = © Vm- 

m=0 

In this case an invariant subspace Vm is generated by vectors: 

and the singular vectors of such a subspace are characterized by the following proposition. 

Proposition 7 Let Vm be an invariant subspace of ft ^ Vx^, with m < mm{Ai, A2}. Then 
a vector v^^ G Vm, 

■m 

V'm ^Y."^ i^'^Z^A, 

is singular if and only if the coefficients Ck satisfy the conditions 

Cfc+i(fc-|- l)(fc - Ai) +Cfc (to- fc)(m - fc- 1 - A2) = 0, V /c = 0,...,TO- 1. (23) 
The proof is based on straightforward calculation using relation (|^): 

m— 1 

= + " ^) + - + fc + 1 - to)} F'^vx, ® F^-^'-^vx,. 

Remark that the conditions (^), satisfied by the coefficients Cfc, coincide with the conditions 
determined in Ref. [|l2j, satisfied by the coefficients of a bilinear differential operator which is 
projective covariant. The system ( p3[ ) was solved in this article and for to < min {Ai, A2} it admits 
a unique solution (up to a constant factor): 

c. = {-lfCt^^^^^^^^, Vfc = 0,...,TO, 
(-Ai)fe 

with Cm the binomial coefficient and {x)i the Pochhammer symbol: 

{x)i = x{x + l)...(a:: + i-l), £ N*, (a::)o = 1. 

By analogy with covariant differential operators, we introduce the bilinear operator Pm defined on 
the subspace Vm of = ® Vx^ by 

Praiv, ®v,)=Y (-1)' Ct ^^~\~^^^^ F'v, ® F-^-'v,, (24) 
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which is analogous to the Gordan operator. 

Conclusion. In the case of = 2 in every invariant subspace Vm with m < min {Ai, A2} there 
is a unique singular vector (up to a constant factor) which is 

= P^{v,, ® V,,) = £ (-1)'= Ci ^"^-^^^^^^ F'^v,, ® F^^-^v,, . 



3.1.2 The case = 3 

Consider now = 3 and = V^i ® Vx^ O with its decomposition on invariant subspaces. In 
this subspace Vm is generated by vectors: 

^2 — 0, . . . , m — fci 

and the singular vectors of this subspace are characterized by the following result. 

Proposition 8 LetVm be an invariant subspace of fl — Va^ (8>Va2 (^Vaj, with m < mm{Ai, A2, A3}. 
Then a vector ofVm, 

m m — ki 

< = E E c,,k, F'^^vx,^F'^^vx,^F^^-''^'''-vx,, (25) 

fci=0 k2=0 

is singular if and only if the coefficients c^^fc^ satisfy the conditions: 

Cfci+l,fc2(fcl + l)(fcl - -^l) + Cfci,fc2 + l(fc2 + l)(fc2 - A2) 

+Cfci,fc2(m -ki- fc2)(m - fci - ^2 - 1 - A3) = 0, (26) 

Vfcl = 0, TO — 1, Vfc2 = 0, TO — 1 — fcl. 



The proof is analogous to that of = 2. As in the case A^ = 2, we note that the system (|2(]|), 
which has to be satisfied by the coefficients of the development (^) in order that be singular, 
coincides with the conditions satisfied by the coefficients of a trilinear differential operator which is 
projective covariant [Q. The system ( p6| ) was also solved in Ref. [|l^ and for to < min{Ai, A2, A3} 
it admits to + I linear independent solutions. In addition, it was shown that the space of covariant 
trilinear operators is generated only by successive applications of covariant bilinear operators. 

Conclusion. In the case of A^ = 3 in every subspace Vm with to < min{Ai, A2, A3} there are 
to + 1 linear independent singular vectors: 

dim (Sing Vm) = to + 1. 



Moreover, using the operator on Q introduced in (27) we can construct a basis of Sing Vm, given by 
the vectors 

{<« = P3,fc (P2,m-fc(wAi ®«A2)«'WA3)}fc=o,...,m- 

3.1.3 General case 

The results obtained for N — 2 and A^ = 3 can be generalized for an arbitrary N. For each N > 2, 

N-l 

consider the space ft = ft^ = fl^^^ CEJ Vajv, where ft^^^ = (g) Va and define on fi^ the operator 

1=1 

PN.mivi ®V2)=y] (-1)' ^"^ , , , FSi ® f "-S2, Wvien''~\v2&Vx,, (27) 



where F = Eilf -P"^'^ is defined by @) on n^-^. Then we have the following. 

Proposition 9 Let Vm be an invariant subspace of il^ with to < mm{Ai, Aat}. Then in Vm 
there is a family of C™^jy_2 linear independent singular vectors, given by the formula 

Vm = Pw.fc ® «A„) , k = 0,...,m, (28) 
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JV-1 

where v^_f. is a singular vector of weight — 2(m — k) in the basis of the subspace SmgVm-k 

i=l 

of n^-\ 

The singularity of these vectors arises from straightforward calculation analogous to the case 
N = 2. 

These vectors are linear independent because w^-fe elements of the basis of the subspace 

Sing Vm~k of f2^^^ and for different k the maximal number of operators F applied on the last 
component is different. 

The fact that by this construction we obtain exactly C'™_|_^_2 linear independent singular vectors, 
can be demonstrated by induction with respect to N: For N = 2 and = 3 this number of singular 
vectors was already obtained. Suppose now that for an arbitrary N the number of linear independent 
singular vectors v^_^ in Vm-k is Then for + 1, the number of linear independent 



singular vectors obtained by the construction (^8|) is 

rn m 

Ey^ni — k \ ^ y^j 

k=0 j=0 

Proposition ^ allows to construct inductively a basis in Sing V™ which has the form 

{PAr,fc„_i(--P3,fc2(P2,fci(wAi ®WA2) 8) fc,=0,...,m ■ (29) 

fcl + fc2 + ■■■ + fcjV-l = "I 

Remark. We point out that in this section we have considered invariant subspaces Vm of fl 
with spin deviation m which does not exceed any of the N weights Xi G N: 

m < min {Ai, Aat} . 

If m is greater than at least one of the weights Ai, then in the set 

{F"^«A, ®...®F""«A„}„,+...+„„.™, 

there are elements for which the number of operators F acting on v\- is greater then A^. These 
vectors are zero because for finite-dimensional irreducible representations of SL{2) we have F^v\ = 
0, Vn > A + 1. Therefore in this case the dimension of the space Vm is less then C™_^j^_-^ and 
depends on the weights A which are less then m. 

Moreover, if m exceeds one or more of the weights A, the dimension of the space Sing Vm is 
also less than C™_|_^_2- For example, in the case iV = 2, it was shown in Ref. that for 

m < rumax = Ai + A2, but m > Ai and m > A2 there are two independent bilinear covariant 
differential operators, to which correspond singular vectors of the type 

^m-(i+Ai) {F^^^^vx, ® vx,) and Pm-(1+X,) (wAi ® F^+^'vx,) , 

which are zero in the considered finite-dimensional irreducible representations of SL{2). For an 
arbitrary N we claim that relation ( p9| ) still gives the basis of Sing Vm but some of the elements of 
the form (E9|) vanish. This fact implies the decreasing of the dimension of Sing Vm- 



3.2 Singular common eigenvectors and Bethe Ansatz 

The basis ( p9| ) of Sing Vm constructed in the previous section is not a basis of common eigenvectors. 
In this section we show that the Bethe eigenvectors constructed in Sec. II using the Bethe Ansatz 
are singular, but their number could be less than the dimension of Sing Vm- More precisely, the 
Bethe equations (p7|), which appear in the Bethe Ansatz as conditions that Bethe type vectors ( p^ ) 
be common eigenvectors of Hi are also conditions that Bethe type vectors be singular. 

Proposition 10 If the distinct parameters wi, ...,Wm are solutions of the Bethe system: 

N m 

E^^+ E ^^-0, Vfc = l,...,m, (30) 

then the corresponding Bethe type vector ^pm{wi, ...,Wm)vo = ^iwi).-.J-{uJm)vo ofVm is a singular 
vector of the tensor product module f2. 
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Proof. Using the commutators, 



[E,T{w)] ^ h{w) and [h{wi), J^{w2)] = 



with h{w) = J2k=i w-z we obtain 



/ 



'Elprn{wi,...,W„i)vo =^ 



k=l k = 1 



\ 



E 

k = l 
k ^ i 



IS a 



Therefore if the parameters wi, ■■■,Wm are solutions of the system (30), then '^^(u'l, 
singular vector. 

The number of common Bethe eigenvectors is equal to the number of distinct solutions of the 
Bethe system (30). If they are C™_|^jy-2 then they form a basis of Sing Vm- For m = 1 and m = 2 we 
can easily show that the number of complex solutions of the Bethe system is — 1 and N{N — 1) /2, 
respectively, but we cannot say how many of them are distinct. For example, in the case m = 1, 
iV = 3, the Bethe equation is equivalent to a second order equation. For some particular values of 
the parameters Xi and Zi, i — 1,2, 3, this equation could have a double solution. 



For arbitrary m and N, the number of solutions of the Bethe system seems to be C, 



m+N-2 



but 



to our knowledge there is no estimate of the number of distinct solutions. In any case, the Bethe 
Ansatz gives a family of singular common eigenvectors which could be completed, if necessary, to a 
basis of Sing Vm with common eigenvectors constructed using the basis (p9|). 



4 The Gaudin model for an arbitrary simple Lie algebra G 

Consider a simple Lie algebra G, of dimension d and rank r. Denote A the root system of G, A_|_ 
the system of positive roots and Aq the system of simple roots. The Cartan basis of G is formed by 
the Cartan generators the generators of positive roots {Ea}^^^ and the generators 

of negative roots {Fa — E^^) ai^i\^- "^^^ commutation relations are 

< a, a > 

[Ea, Ep\ = Na,pEa+i3, V Of, ^5 € A, such that a + /3 e A. 

The Killing form of G defines a symmetrical, G-invariant, bilinear form < , > on G, which allows 
us to identify the Cartan sub-algebra /C of G with its dual, by the isomorphism a{H) =< iJ^, H >, 
for all iJ G /C. The scalar product on K, induces a scalar product on its dual, which is nondegenerate, 

< a,P>^< Ha, Hp >, 

normalized such that < ai,ai >— 2 for a long root a/. 

All generators of G are orthogonal with respect to the bilinear form < , > except 

<H„H, > = ^x, Vz=l,...,r, 
2 

<Ea,Fa> = , VaGA+, 

< a.a > 

where x — Y^l=i '^fi^i) the square of the length of a long root. Then, for the Cartan basis we 
can define the dual basis with respect to this bilinear form: 

H-ln- g, _ <a,a> ^ ^ - _ <a,a> ^ 

X 2 2 

such that < Hi, Hj >— Sij and < Ea, Ep >—< Fa,Fp >— Sap- 
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For this algebra, consider N finite-dimensional highest weight G- modules, V\-^,..., V\j^, with 
dominant integral highest weights Ai,..., Xn and highest weight vectors waivj ^An- The tensor 
product module of these G-modules is the space of physical states: 

n^Vx, ® ...(8)yA„. 

Recall that the Lie algebra action on this tensor module is defined by (|l|). 

The fundamental building blocks of the generalized Gaudin Hamiltonians are the quadratic 
generalized Casimir operators, defined on Q as 



= E^a^^^i^'^ Vz,j = l,...,iV, (31) 



a=l 



where {Ia}a=i d^^^ basis of G and \ Ia\ its dual with respect to the scalar product, < Ia,Ib >= 

' L J a—l,d 

Sba- We recall two fundamental properties of these operators, which will be useful in this section. 

(1) w''^-* are independent of the choice of the basis in G. In particular, in the Cartan basis they 
take the form 

^ 1=1 QeA+ 

(2) Lu^'^^'f commute with the action of the algebra on the tensor module: 



= 0, VX e G. 



The operators w'*-'-' do not commute between themselves, but we can construct N linear combi- 
nations which commute. 

Lemma 11 The operators uji : 11 — > fi, i — ■■■■,N, 

N 



E 



commute: [uJijUjj] ~ 0, Vi,j = 1,...,N if and only if the coefficients Cij satisfy the equations 

CikCjk - CikCji - CjkCij = 0, yiy^j^ky^i,i,j,k^l,...,N. (32) 

In particular, if zi, zat are distinct complex parameters, the coefficients Cij = j't^ fulfill the 
conditions (p2[). With these coefficients we can construct the following N operators on (l: 



N 

H^= y -^uj^'^\ (33) 



which are the generalized Gaudin Hamiltonians. They preserve the properties (1) and (2) of the 
operators u!^^^\ More precisely, if we consider the Cartan basis in G, we have 

[^^,Ha] = 0, (34) 
[7^„EJ = 0, Vi = l,...,iV; Va e A+, (35) 
[H„F„] = 0, (36) 

where Hq — J2iLiHa\ = J^iLi^aK — J^iLi^a'' are operators on the tensor product 
module, defined by (|l|). The generalized Gaudin Hamiltonians commute between themselves, but 
only A'^ — 1 of them are independent. Due to the property ( ^ ) we can complete this system by 
other r operators, the generators Hq of the Cartan sub-algebra, which commute with all Tii. For 
this family of independent and commuting Hamiltonians there is a complete system of common 
eigenvectors in fl. 
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4.1 The structure of il 

Since the space is a tensor product of N highest weight representations Vx of G, it is useful to 
recall some basic results concerning the theory of highest weight representations of Lie algebras. 

Such a representation is completely determined by a highest weight vector v\, on which the 
action of the algebra is given by 

HaVx =< a, A > Vx, EaVx = 0. 

As in the SL{2) case, the representation space Vx is linearly generated only by monomials in gen- 
erators of negative roots: 

{v2\::2 = Kl-Kl vx}o..e^^. n.en, (37) 

but in this case, the monomials are ordered with respect to the roots, as showed by the Poincare- 
Birkhoff-Witt theorem. These vectors are eigenvectors for the Cartan generators: 

Ho.vl\::Z = <a,A>-^n^<a,/3> (38) 
V /3eA+ / 

and the action of a positive root generator on monomials w^" is 

Jga<° = f 2 ^ " ^ - n„ + vl--\ (39) 
\ < a, a > / 

Concerning the dimension of a highest weight representation of G, recall that a weight A is a 
dominant integral if 

r„ ^ 2 ^ ^' " ^ £ N, VaeAo. 

< a,a > 

If the highest weight A is not dominant integral, then the representation is infinite dimensional and 
irreducible. If the highest weight A is dominant integral then there is an invariant subspace,[T^ 
generated by vectors w^J^ ' ^j, with tiq > + 1. The quotient representation, denoted Vx is an 
irreducible, finite-dimensional representation of G, generated by ordered monomials: 

f ni...ni p^l Zp^i \ 

\'^ai...ai ~ ^ai •■■-'^^ai ^\ } a, eA+ , nj=0,. . . ,ra^ ■ 

As in the SL{2) case, the space is a finite-dimensional tensor product module of G, completely 
determined by the vacuum vector vq = vx^ (8) ... iX> ^nd generated by ordered monomials: 

with 1 < ki < N. Note that for finite-dimensional representations, m can vary between and 
a maximal value TO„ax = X]i3=i X^aeA^ ^a- Note also that the action of the Hamiltonians on the 
vacuum vector is 



N 



' . r-',. z, - Zj I 

As in the case of the Lie algebra SL{2), the Cartan operators have the particular role to 
give a gradation of the tensor product module: 
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N 



on weight subspaces V^, of weight /i = Ai — jj, generated by eigenvectors: 



N 



-.1 j=i 



'vq, 



with 7i, ...,7m € A_|_ and 1 < fc^ < A^. Since a positive root can be written as a sum of simple 
roots, with positive integer coefficients, in a unique way, the weight subspaces wih be labeled by 
a family of simple roots: V^i = Vq.j...q.^, with ai...as G Aq, not necessarily distinct and satisfying 
Ej" 1 7j = Efc=i "fc- Hence 

^ ^ ^Ctl ■■■Cts • 

From the explicit form of the generalized Gaudin Hamiltonians ( p3| ) and from the action (|37|)- 
( |39| ) of the Lie algebra on each module , it follows that the weight subspaces Vai...as ^"^^ invariant 
under the action of any Hi. Therefore, in each subspace Va^...a^, we can construct a basis of common 
eigenvectors of Hi- 

Unlike the case of the Lie algebra SL{2), there is no result concerning the dimension of the 
invariant subspaces Vai...a^ and neither for Sing Vqi...q,^. However, using the properties ( |35| ) and 
(|36|), the recursive procedure presented in Sec. HI can be generalized to construct in each subspace 
Vai...as 8- family of common nonsingular eigenvectors, by applying operators Fq. on vectors of the 
sub-spaces Vq,j...q._jq._|_j...q,^. Concerning the subspace of singular vectors in Vai...as do not have 
a generalization of the result presented in Sec. HI A, for the construction of a basis of singular 
vectors. The only result Q concerning this problem is a generalization of the Bethe Ansatz, which 
we discuss hereafter. 



4.2 Common eigenvectors in Va 

Consider a simple root a. The subspace Vq of weight X)t=i is generated by vectors I Fa^^ \ 

^ J fc— l,...,A^ 

A Bethe vector in this subspace is defined as an expansion on the basis, with rational coefficients 
depending on one complex parameter w: 

V'i(w,a)wo = — - — Vii = T{w,a)vQ. (40) 

2—1 

To give the action of a Hamiltonian Hi on this vector it is useful to calculate the commutator: 
[Ht,T{w,a)] = T{'w,a) — — V — — 

W — Zi W — Zi \ ^ — ' W — Zj 

„ 2 w — Zi w — Zi w — Zj I 

Applied on the vacuum vector all the terms of the last sum vanish and the action of a Gaudin 
Hamiltonian Hi on the vector (|4^) is 

n-i I I \ I < Aj, Aj > < a, Ai > \ (^<a,Xk>\ F^^ 

Hence we have the following. 

Lemma 12 Given N complex numbers {zi}^^^ ^ and dominant integral highest weights {Xi}^^^ ^, 
the Bethe vector tpi(w, a)vo is an eigenvector for all Gaudin Hamiltonians: 

Hi ipi{w, a)vo = sl{w, a) ^i{w, a)vo, V i 1, N, 
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if the complex parameter w satisfies the Bethe equation: 



N 



^ < g, Afc > ^ p 



fc=i 



W- Zk 



(43) 



The eigenvalue s]{w,a) of the Hamiltonian Hi depends on the solution of this equation: 

N 



< Xi, Xj > < a, Xi > 



i = 1,...,N. 



4.3 Common eigenvectors in Vai,a2 

Consider two simple roots ai and a2, which are not necessarily distinct. The subspace Vq.j.q2, 

of weight y^^i Xi — ai — a2 is generated by vectors \ F^,^'' f!^^^\^\ , with 

two generators of negative roots applied on two distinct components of wq, but also by vectors 
{(i^QjF„2)('^)wo}j,^i ^ and {{Fa^Fa^Y^^'vo} ^ with two generators of negative roots applied 
on the same component of wq. A Bethe vector in Vax,a2 is defined as an expansion on all these 
vectors, with some particular coefficients depending on two complex parameters wi,W2- 



ip2(wi,ai;w2,a2)vo 



N N 

E E 

fci = l k27ik 



p(fcl) p(fc2) 



1 ^ 

-^E 



(F F 

OLl^ 0.2} 



1 ^ 

-^E 



(F F "l^^) 



-Vo- 



As element of the representation space V\^. , one of the two last terms is not a good ordered monomial 
and we must write, for instance, (Fq^Fq-j )'^'^) = (Fq^Fq^)''^) + [F^^ , Fq J Then this vector can be 
also written in the following form, using the operators 



■)p2{wi,ai;w2,a2)va = T {wi, ai) T {w2, a2) vq + 



Wi — W2 



-T{wi, [Fa^,FaJ)i;o- 



(44) 



In order to compute the action of a Hamiltonian TLi on the vector ^/;2, we use the result ( p] ) and 
|) to obtain 



Uiii)2{wi,auW2,a2)vo = s'^{wi,W2) ■4'2iwi,ai;w2,a2)vo 



F^} 1 

— J^{w2,a2) H 

Wi — Zi Wi — W2 Wi — Z. 



a-L ; ^ a2 



-f2 T{wi,ai)- 



F*- 



^ ai , ^ a2 



W2 - Zi Wi - W2 Wi - Z, 



«0 



«0, 



where 



<A.i,Aj> <ai,Xi> <a2,Xi> 



Wi - Zi 



N 

E 

N , . 

E< ai, Afe > < ai, a2 > , . ■r^ < q;2, > < ai, 012 > 
,„ ^„ + ,„ and /2 = 2^' 



W2 - Zi 

N 



fc=l 



Wl - Zk 



W2 — Wl 



k=l 



W2 - Zk 



Wl — W2 



Note that the action of Tii on the second term of the vector (Q) must be calculated separately, 
because ( |4^ ) does not hold for the positive root ai + a2 which is no more simple. It follows from 
these considerations the following lemma. 
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Lemma 13 The Bethe vector ip2{wi,ai;w2,a2)vo is an eigenvector for all Gaudin Hamiltonians, 

'Hiip2iwi,W2)vo = Si V'2(wi,ai;w2, 02)^0, 

if the parameters wi , W2 satisfy the Bethe equations: 

N , N , 
/i = > \ =0, /2 = ^ = 0- 45) 

4.4 Common eigenvectors in ... 

For m simple roots ai, not necessarily distinct, the subspace Vqj^..._q,^ of weight X^ili ~ 
is generated by ordered monomials F^'^''\..Fj^"''vq with ki, fc„ = 1, N and 71, ...,7n 
positive roots with J2"=i 7* = SfcLi Q^fe- I^i this general case it is difficult to define an appropriate 
Bethe type vector, with rational coefficients on the basis of Vai,...,a„^- In Ref. Q| a recursive 
procedure was proposed to define such a vector: 



)vo = ilJra~i{wi,ai; W^n-l, am-i)^ {Wm, am) Vq 

m — 1 ^ 

+ i'rn-i (wi,ai; ...■,Wj, [Fa.,Fa] ; Wm_i,a„_i) vq, (46) 



with ipi^w, a) = ^ {w, a). However, to our knowledge there is no proof that the action of a Hamil- 
tonian TLi on the vector ipmVo is 

m 

Ti-i '4>miwi,ai; Wm, am)vo = s"^ -(/'m(wi, ai; Wm, Q;m)'yo - ^ /rV'm'^o, (47) 

fe=l 

with ■0^ some vectors in Vai,...,am ^-^d 

- 1 ■ / - Zi Zj W}^ Zi 

TV . m 

Such a recursive construction of generalized Bethe vectors seems to be appropriate for inductive 
calculations, but for the action (|47| ) of Hamiltonians on these vectors such a calculation raises some 
problems, as explained in Ref. ]10( |. Similar problems occur if we intend to prove inductively that all 
Bethe eigenvectors are singular. For a small number of roots this can be done by direct calculation. 

We call ti* G a singular vector of if all the generators of simple roots Eq act trivially on v" : 

E^w^ = 0. 

Note first that [Ef^, Fa] = Sa/i ^aa> simple roots a and (3. 

Consider now the weight subspace Va- If /? is a simple root we have 

2 ^ Hi^^ 



< a,a > ^-^ w — Zi 

Hence, 

E^Vi(w,q;)wo = (5q/3 ( ^ ]vo, V/3 e Aq. 

< a,a > w — Zi j 

Therefore, if w is a solution of the Bethe equation (|4^), then i^i{w^ a) wq is a singular vector of Va- 



N , 

< a, A; > 
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In order to calculate the action of a generator of simple root on a generalized Bethe vector 
(H) of Vai,Q2, we use (HH) but also 

<C Oil Ct2 ^ 

[E,f3,^pi{w,[Fai,Fa2])] = 2— — {5/3a2-^(lf,ai) - (5/3qi ^(w, 02)} , 



<f3,f3> 



to obtain 



E;3V'2(wi,ai; W2,a2)wo = S|3a^ 



2 j < Aj,ai > ^ < Q;i,a2 > 

< Qfl, Q!i > 1 Wi — Zi W2 — Wi 



N 



E<Xi,a2> <ai,a2>. . , 
~. : ^ ~. ~. — r ^{wi,ai)vo. 



<a2,a2> \ ^ W2 - Zi Wi - W2 



Therefore, if {wi,W2) is a solution of the Bethe equations (p^, then '02(^1, ai; W2, 0^2) is a singular 
vector of Vq^.q^ . 

For Bethe vectors (^|) depending on more than 3 simple roots we claim that the commutator of 
with ijjjn is given by 



2 J ^ i^W) 

\^l3,i'ra{wi,ai;...;wm,a.„^)\ = y^'^/Sa; < V^m-i(---;5^7ai; --Oy^ ■ 

^ < a,:, a,; > ' ^ 



E 



' 1pm-l{---Wi,ai...) 



where ■>jjm-ii---Wi^i---) denotes '0m-i(---'f«i-i: cti-i; Wi+i, cti+i---)- As for the action ( ^ ) of the 
Hamiltonians, there is no inductive proof of this relation. The action of on the Bethe vector is 
then 



E/3V'm(u'i, ai; 'Wm,am)vo 



2S 



N , m 



I /-^ 1,, , _ V : ^-^ 



fc = 1 



V'm-l(--.W^7ai---)^0- 



Hence, the Bethe vector is singular if the parameters w satisfy the Bethe equations /™ = 0, V fc 



1, ...,m. 



5 Conclusions 

In this article we recall the Gaudin model associated to the Lie algebra SL{2) and its partial 
diagonalization by the Bethe Ansatz. 

We give a general recursive method to construct a basis of common eigenvectors in each invariant 
subspace Vm- Knowing a basis of Vm-i, we construct a family of nonsingular independent common 
eigenvectors in Vm- This family is completed to a basis of V™, by a basis of the subspace of singular 
vectors of Vm- 

In order to describe the subspace Sing Vm, we establish a relation between singular vectors and 
covariant differential operators. This allows us to construct a basis of Sing Vm using an analog of 
the Gordan operator. On the other hand we show that the Bethe Ansatz gives a family of singular 
common eigenvectors. If Bethe equations have a maximal number of distinct solutions, then Bethe 
eigenvectors form a basis of common eigenvectors in Sing Vm- 

We discuss also the generalization of this method to the case of an arbitrary Lie algebra. We 
recall the generalization of the Gaudin Hamiltonians and of the Bethe Ansatz. The generalized 
Bethe vectors are defined recursively, but this definition allows neither to calculate the action of 
the Hamiltonians, nor to prove that Bethe eigenvectors are singular. For a small number of simple 
roots we prove by direct calculation that the Bethe equations are conditions that Bethe type vectors 
be singular and common eigenvectors. The recursive method to construct nonsingular common 
eigenvectors could also be generalized to the case of an arbitrary Lie algebra. 
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